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We consider A-type model of the Bose-Einstein condensate of sodium atoms interacting with 
the light. Coefficients of the Kerr-nonlinearity in the condensate can achieve large and negative 
values providing the possibility for effective control of group velocity and dispersion of the probe 
pulse. We find a regime when the observation of the "slow" and "fast" light propagating without 
absorption becomes achievable due to strong nonlinearity. An effective two-level quantum model of 
the system is derived and studied based on the su(2) polynomial deformation approach. We propose 
an efficient way for generation of subpoissonian fields in the Bose-Einstein condensate at time-scales 
much shorter than the characteristic decay time in the system. We show that the quantum properties 
of the probe pulse can be controlled in BEC by the classical coupling field. 
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I. INTRODUCTION 



The problem of atom-field interaction represent one of 
the major areas in the modern physics and quantum op- 
tics, particularly [l], |2|. A special interest is regarded 
to effects of propagation of light inside highly nonlinear 
atomic media where a large nonlinearity is achieved with 
appropriate choice of the form of the atom-field interac- 
tion 0, 13 ■ Some early works [a, Q on interaction of a 
resonance atomic system and laser field demonstrated a 
possibility to achieve significant difference between phase 
and group velocity of the light pulses. In the experiments, 
Lctokhov and Basov showed that while the front-edge of 
the pulse generated the inversion in the system of two- 
level particles resulted in the sloping of the front-edge, 
the back process of re-emission of the absorbed energy 
gave rise to some steeping of the back-edge of the pulse. 
Hence, the pulse shape deformation resulted in signifi- 
cant delay of its registration on exit from the resonance 
media that was observed experimentally. The works mo- 
tivated intensive study of so-called the " slow" and " fast" 
light. 

An important step in the development of the theory 
was the theoretical prediction |j| and experimental ob- 
servation |3 of positive and negative delays of picosec- 
ond pulses propagating without pulse shape deformation 
inside a crystal. But, very high level of optical losses was 
limiting the effect to a relatively small magnitude. Ba- 
sically, there are several different ways to overcome the 
difficulty The first idea exploits the shortening of the 
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pulse duration to lengths, which would be much smaller 
than the relaxation times in the medium that provides 
necessary conditions for a generation of optical solitons 
[9j. Another approach is based on a three-level A-scheme. 
In that case one of two laser beams is a strong coupling 
field developing a transparency window in the medium 
and the second probe field propagates through the res- 
onant system without absorption and with unchanged 
pulse shape |10| . 

A sketch of the energy levels in the A-scheme for a sin- 
gle sodium atom is shown in Fig. ^ A classical beam 
with the central frequency u c couples the levels |1) and 
|3), and the probe pulse with central frequency uj p cou- 
ples 1 3) and |2) levels such that the three levels form a 
A-type configuration. If the intensities of the classical 
and probe fields are comparable and the relaxation rates 
of the levels |1) and |2) are negligible then a state of 
an atom may be considered as quantum superposition of 
two lower levels |1) and |2). It corresponds to the effect of 
so-called coherent trapping of lower states [ill Il2l | . An- 
other type of dynamics may be realized if the intensity 
of the probe pulse is much smaller than the intensity of 
the classical field and the atoms initially populate the 
lower level |2). The regime is called electromagnetically 
induced transparency (EIT). In the paper we concentrate 
on the later. 

Systems in the EIT regime are extensively studied in 
the literature. Some interesting results are referred to: 
the studies of the dispersion of the atomic medium in 
a linear response regime [l3j, semiclassical estimates of 
the time delay of the probe field inside atomic cells in a 
thermal equilibrium [l4[ , EIT in doped crystals [15| , non- 
linear optical parametric processes in resonant double-A 
system [l|j. Wang and colleagues carried out an impor- 
tant experimental measurements of the Kerr nonlinear 
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FIG. 1: The energy level A-scheme of 23 Na atoms. 



index of refraction in a three-level Rb atomic A-system 
using an optical ring-cavity 17]. They found that the 
Kerr nonlinearity can reach very large and negative val- 
ues. 

One of the dominating problems in EIT was the 
spectrum broadening regarded to thermal effects. To 
overcome the difficulty, a highly coherent atomic Bose- 
Einstein condensate (BEC) is used [18|,[l|j. Hau and co- 
workers calculated linear dispersive properties of the sys- 
tem and demonstrated a possibility to observe the " slow" 
light effect in the ultra-cold vapor of sodium atoms loaded 
into a magnetic trap at the temperature of transition 
to the condensed state 20]. Based on their results we 
consider EIT effects in the Bose-Einstein condensate of 
sodium atoms. Two hyperfine sub-levels of sodium state 
3 2 <Si/2 with F = l,F = 2 are associated with levels 
|2) and |1) of the A-scheme, correspondingly. An ex- 
cited state 1 3) corresponds to the hyperfine sub- level of 
the term 3 2 P 3/2 with F = 2 (see Fig. [TJ :21J. The en- 
ergy splitting between the levels |1) and |2) is denoted 
by cj 12 /(2tt) = 1772MHz, the transition |3) -> |2) cor- 
responds to the optical frequency oj/(2tt) = 5.1 • 10 14 Hz 
(A = 589nm) l^. 

The condensate is considered being placed inside a 
confocal resonator. A strong coupling beam propagates 
along the optical axis of the resonator maintaining the 
transparency window in the condensate due to the tran- 
sitions |1) — > 1 3). According to the conventional EIT 
regime the atoms mainly populate the lowest level |2). 
We assume that the resonator mode has an eigenfre- 
quency u> c . The coupling laser beam develops large po- 
larization in the condensate resulting in significant non- 
linear susceptibility of the medium on the frequency u> p . 
A probe pulse with the appropriately chosen polarization 
direction propagates transversely to the resonator optical 
axis. It obeys high nonlinear interaction with the conden- 
sate. In order to study a competition between linear and 
nonlinear processes regarded to the transition |3) — > |2) 
induced by the probe pulse, the linear and third order 
Kerr-like terms in the expansion of atomic susceptibility 
are considered. An effective control on the absorption 
of the probe pulse and its group velocity is realized by 
detuning S of the coupling field frequency from the reso- 



nance. 

In the following section we find conditions when the 
group velocity of the probe pulse becomes very small or 
very large. We also note a region where the absorption 
in the condensate is almost zero what we have called as 
nonlinear compensation of optical losses. In the third 
section, the interaction between sodium BEC and the 
probe pulse is described by an effective quantum Hamil- 
tonian. We explain the way to reformulate linear three 
level A-scheme in terms of nonlinear quantum model of a 
two-level particle interacting with quantized field and to 
solve the model. We apply the su(2) polynomial defor- 
mation formalism in order to develop the perturbation 
theory. In the next section the perturbation theory is 
considered up to second order. The last fifth section con- 
cerns with some nonclassical effects in quantum statis- 
tics of the photons in the probe pulse described by the 
effective Hamiltonian. In the conclusion we discuss some 
experimental perspectives of our results. 



II. NONLINEAR COMPENSATION IN SODIUM 

BEC. 

Applying the formalism of slowly varying field ampli- 
tudes |23J in the rotating wave approximation we write 
the Hamiltonian describing the interaction of three-level 
atoms with two laser fields in the following form |24|: 

H a =lu 12 \1)(1\+lu\3)(3\ 
- 5l |3)(l|e-^(*-t)- ffl |l)(3|e^(*-t) 
-g 2 \3)(2\e- iu ^ t -^ - g* 2 \2) {^e^^- ^ . 



(1) 



For simplicity, we let ft = 1 in the paper. In Eq. QJ the 
coefficients g\^ 2 determine single-photon Rabi frequencies 
and are defined as follows 
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|M3l|^c, 32 = \/J>32\Ap. 



(2) 



Here, /Xjj is the atomic dipole momentum, A c ^ are the 
slowly varying coupling and probe field amplitudes, cor- 
respondingly. 

We denote single-atom density matrix by p = 
J2ti 7>=i Pij\i)(j\- The time evolution of the density ma- 
trix is described by the Liouville equation of motion |25j 



di 



- 7l2 (|l}(l|p-2|2)(l|p|l)(2| + p|l)(l|) 
-732(|3}(3|p-2|2)(3|p|3)(2| + p|3)(3|) 

-73i(|3}(3|p-2|l)(3M3)(l|+p|3)(3|). 



(3) 



Here, the constants 7y determine the rate of spontaneous 
decay from the levels \i) to \j) in the A-scheme. In gen- 
eral, to consider the space-time dynamics of the system 
of fields and atoms interacting in the resonator we have 
to add to Eq. (|3J) the Maxwell equations 
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Here, c is the velocity of light in the vacuum, E is the to- 
tal amplitude of corresponding fields. The vector P is the 
polarization of the condensate induced by the coupling or 
probe field, correspondingly. 

In the adiabatic limit, when the variation of the Rabi 
frequency <?i is very small (26| a self-consistent problem 
of Eqs. © and (|3fa[) may be reduced to a less complicated 
one. In that case the system of equationscan be solved 
separately for the medium and the fields. Averaged the 
density matrix elements over the rapidly oscillating phase 
of the fields we represent p in the following form 
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A total polarization vector of the BEC coupled to the 
probe electromagnetic field is given in the form P — 
p(0 + p(m) [27^ Here p(J) = eo £(i) . g describes the 
linear contribution, and 
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together with the relation pij = p*^ A bar in the matrix 
elements p denotes the averaged quantities. Substituting 
definitions Eq. Q) into Eq. (JSJ) one obtains equations of 
motion for p matrix elements 



is the nonlinear response on the external fields. The hat 
denotes tensors. Using well-known relation for the po- 
larization induced in a resonance medium P = =^-^32/932 
|2| and Eq. J7J we find first and third order nonlinear 
susceptibilities of the Bose-Einstein condensate: 
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Here, A = uj — uj p and 5 = u> — W12 — ^ c - 

Since we are interested in nonlinear interaction of the 
BEC with the probe pulse in the dipole approximation 
we only consider the explicit dependence of P32 on the 
Rabi frequency 32 of the probe field, 



Here Af is a number of atoms in BEC, V is the volume, 
(5)and Af/V is the density of the condensate. Earlier, a 
similar form for the linear susceptibility \ °f BEC in 
the EIT regime was obtained in the paper 20]. The 
nonlinear part x of Kerr type was studied in the regime 
of giant nonlinearities induced in a cyclic process of A- 
type interaction between optical fields and "hot" 87 Rb 
atoms in atomic cells inside an optical ring cavity |17| . 

The permittivity of Bose gas corresponding to the 
probe field and including linear and nonlinear terms, 
reads IH 



P32 
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= 1 + X (1) + 7X (3 W 
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where p 32 = denotes initial polarization in the con- 
densate, which is zero for the sodium. The coefficient 
P32 corresponds to the stationary solution of the system 
Eq. (jSJ in the linear approximation and it is responsible 
for the linear susceptibility of the medium. The nonlin- 
ear corrections p 32 and p 32 determine resonant nonlinear 

atomic susceptibility. The calculation shows that p 32 it 
is negligible comparing with the Kerr type nonlinearity 

p 32 \ Hence, we neglect the second order correction p 32 
in the paper. 

We study the dynamics of |2) and |3) levels of sodium 
atoms interacting with the probe field. We characterize 
the transition |2) — > |3) in terms of the effective coupling 
constant associated with the dipole matrix element p^- 
Considering all atoms being initially in the state |2), i.e. 
P22 = 1, Pn = P33 = we find from Eq. JSJ 



Hence, using the relation e p = I n p + i^- ) and Eq. i|ll|) 

we find the refraction index n p and the absorption coeffi- 
cient rj p in first order with respect to the intensity of the 
probe field 
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In the work, the Bose-Einstein condensate is described 
by the A-scheme in near to resonance regime. The con- 
centration of sodium atoms in the condensate Af/V = 
3.3 • 10 12 cm~ 3 is taken from [2(j and the coupling field 
intensity I c — 552^. Taking the dipole matrix element 
|Msa| =22-10 """ 
definition A c = 
stant <7i/(27r) = 21.4MHz. The probe pulse is assumed 



30 C • m from |28| and making use of the 
2"H we calculate the coupling con- 
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FIG. 2: The refraction index n p and the absorption coeffi- 
cient rjp of 23 Na BEC as functions of the detuning ■£-. 



FIG. 3: The group refraction index n g as a function of the 
detuning ^-. 



to have a time-length of 1/us, the laser waist inside the 
BEC is d = 3.7/im. The intensity of the probe pulse 
Ip = 80-^— j- corresponds to 25 photons in average. 

Due to atomic coherence in BEC, in absence of the 
Doppler broadening the decay rates 731 and 732 of the 
level 1 3) can be estimated by natural (spontaneous) width 
of transitions from the upper levels of sodium atoms. 
Taken from the paper J2g] the lifetime on the upper level 
is T re i — 16.3ns, so we assume 732/277 = 731/271 = 5MHz. 
The decay rates of transitions between the hyperfine lev- 
els |1) and 1 2) is 712/(277) = 38KHz according to @. 

In Fig. |3 we plot a typical frequency dependence of the 
nonlinear optical refraction index and absorption coeffi- 
cient as functions of detuning A. Giant nonlinear refrac- 
tion index formed in the condensate with appropriately 
chosen detuning of the probe pulse is exploited in the 
paper to demonstrate the possibility of generating sub- 
poissonian statistics of photons on very short time-scales. 
It is also worth noticing that the realization of a negative 
rip ' is similar to the effects observed in [lT], and it may 
have some important physical applications. 

On the other hand, the alternation of regions with pos- 
itive and negative absorption coefficient corresponds to 
the regimes of effective nonlinear reduction or amplifi- 
cation of the probe field intensity in the BEC. In both 
cases, the energy is transformed between the coupling 
field and the probe pulse. The case of zero absorption 
j] p = may be characterized as a nonlinear EIT. Notice 
that the regime of complete absence of optical losses in 
the system falls at the region of very large values of n p . It 
opens even more extensive perspectives to generate non- 
classical light in the system. In generally, varying the 
intensity of the coupling and probe light, and the detun- 
ing A we can control the parameter T in Eq. (Ill|) and 
change linear and nonlinear relative impact on the probe 
pulse in the condensate. 

To understand the effect of " slow" light in the system 
we study a frequency dependence of some characteristics 



of the probe pulse envelope, 
group refraction index n g 



In Fig. we plotted the 
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A point A = 
of the exact resonance corresponds to the EIT regime 
[lfl l characterized by low losses and large n g . Under cho- 
sen resonance conditions, the group velocity v g = c/n g 
is 2000to/s while the losses are r\ p = 242m, -1 . So, the 
regime corresponds to the effect of "slow" light. An in- 
teresting point in Fig. [5] is a zero of the group refraction 
index n g that indicates an uncertainty of the group ve- 
locity v g of the probe pulse. The uncertainty depicts a 
possibility to observe super luminal velocities. This ef- 
fect together with the negative group refraction index 
is explained in the literature [3, l29| by the presence of 
plane waves of different frequencies, which appear in the 
medium long before the pulse peak enters into it. When 
the group velocity is negative (n g < 0) the pulse gen- 
erates two anti- propagating pulses in the BEC. One of 
them propagates backward (negative velocity) and pre- 
vents the peak from traveling forward in the medium. 

In the following sections we concentrate on some ef- 
fects in dynamics of the weak probe pulse. These effects 
are induced in the system of highly correlated three-level 
particles by relatively strong classical coupling field. We 
show how the later can control some quantum properties 
of the probe pulse. 



III. THE EFFECTIVE TWO-LEVEL QUANTUM 
MODEL 



We are interested in quantum properties of the probe 
field and define the field in terms of creation and anni- 
hilation operators a\a. Considering only two levels |2) 
and 1 3) of interest for the dynamics of the probe field we 
work in the EIT limit, i.e. most of the atoms is concen- 
trated on the lower level |2). We assume that the clas- 
sical field is strong and describe its influence in terms of 
effective coupling constant for the probe field. The effec- 
tive two-level Hamiltonian in dipole and rotating-waves 



approximations reads 

H = Lu p (a^a + S 3 + ^-) + AS 3 + k 1 (aS+ + a*S-) 

+k 2 (a^aa 1f S- + aa 1f aS+). (15) 

The operators S± , S3 describe total dipole momentum 
corresponding to the transitions |3) «-» |2) for the atoms 
in BEC. The first and the second terms give the free 
energy of the probe field and atoms. The third term de- 
scribes linear contribution into the interaction between 
the field and two-level particles in dipole approximation. 
It has the typical form of so called the Tavis-Cummings 
|30J or Dicke model ^. The last term in Eq. iJTSJ de- 
scribes nonlinear processes due to the presence of strong 
classical coupling field, and it depends on the intensity 
of the probe pulse. Both coupling constants are defined 
below. 

An interaction between the BEC and the probe field 
is determined by the effective matrix element of atomic 
dipole momentum between |2) and |3) levels. Due to 
induced nonlinearity the matrix element depends on the 
intensity of the probe field described by the operator a' a 
and on the intensity of the coupling field given by a c- 
number parameter. We define overall coupling constant 
in the dipole approximation as follows 
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(16) 



The generators Mo , M± of the algebra are realized as 
follows 

M = aia ~ S \ M+ = (1 + ^a)a^S (19) 

2 k\ 

and M_ = (M+y . These generators Mq,M± satisfy ba- 
sic commutation relation for any PAE 



[M ,M ± ]=±M ± , 
and commute with the operators 



(20) 



M = a^a + S 3 + r, S 2 = Sf + - {S+S- + SS+) . (21) 



Hereafter we use the same notation M both for the 
Casimir operator and its eigenvalue, if no confusion 
arises. It is known that the eigenvalues M,r(r + 1) of 
the operators in Eq. (|21|l parameterize the PAE in ques- 
tion. We thus denote this PAE by M.M, r - The structure 
polynomial of W^M,r can be expressed in the form 
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(22) 



The parameters of this structure polynomial are 
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where p 32 is given in Eq. @), and fco = /Z32 

the single-photon Rabi- frequency in the Dicke model |32|. 
The form- factor in Eq. (|16|l satisfies the condition that at 
zero coupling field it approaches unity. This behavior of 
the form- factor results from the fact that when g± s=s the 
induced nonlinearity for the probe field should vanish. 

From the definition Ea. l|l()l) and the expansion of p 32 
in Eqs. iJBJ, J7|) we find the coupling constants k\ and k 2 



k x = k L t , k 2 






(17) 



The parameters L;( n ;) denote linear and nonlinear form- 
factors (compare with 33]), correspondingly 
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As far as the phase of the probe field is arbitrary we can 
choose it in such a way that k\ would be purely real with 
Li = \Li\. In that case the nonlinear coupling constant 
would be defined with L n i = e _iarg(Li) L„i. 

According to this formulation we are able to dynami- 
cally control the rates of the transitions induced by the 
quantum probe field A p with help of the classical field 
A c . Changing the latter one can reach qualitatively dif- 
ferent regimes of the quantum dynamics described by the 
effective Hamiltonian Eq. I|15|) . 

The effective Hamiltonian can be expressed in terms of 
fifth-order polynomial algebra of excitations (PAE) [34| • 
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Notice that the root go is a degenerate root of order 2. 

Following standard procedure |3^| we describe physi- 
cally interesting finite dimensional irreducible represen- 
tation (irrep) of Mm,?-- In our model the parameter r has 
the meaning of collective Dicke index (an analog of the 
orbital momentum) of the system of Af two-level parti- 

cles. This index runs from e(Af) = — ^j-^ — to =%■ with 
unit steps, while M can be any natural number including 
zero that follows from the definition Ea. (j21|> . Physically 
interesting irrep of M.M,r is characterized by two roots of 
the structure polynomial where it takes positive values 
in between. The roots 91,2,3 do not depend on the ra- 
tio y~ and are ordered according to the relation between 
M and r. For M > 2r q x < q 2 < q 3 , and for M < 2r 
q 2 < gi < q 3 . Since the root q^ is of order 2, its position 
on the real axis doesn't influence the region where the 
structure polynomial is nonnegative. Hence, if M > 2r, 
the irrep is called a remote zone and ps(x) is nonnega- 
tive in the interval [q 2 , q 3 ] , whereas if M < 2r the irrep 
is called a nearby zone and ps(x) is nonnegative in the 
interval [91,(73]. Notice that the region M ^> 2r is usu- 
ally called the strong-field limit and the region M -C 2r 
is usually called the weak- field limit. The case M = 2r is 
of special kind and the corresponding irrep is called the 
boundary zone. 



Depending on the value of the ratio j^ we might have 
two different situations. If qo doesn't belong to the in- 
terval where ps(x) > 0, the polynomial is approximated 
by the parabolic function relatively well that correspond 
to the method described in our previous paper [3-'l| . In 
case go belongs to the interval of positiveness of ps(x) 
we must introduce some changes to the approach. But 
for the physical system considered here the ratio takes 
large negative values and the total number of excitations 
M may be considered being smaller than the ratio, i.e. 
M < 1 1 1 1 — 1. Therefore, q is always larger than q 3 and 
we can use the algebraic approach developed for conven- 
tional Tavis-Cummings model. 

In the physical situation being studied here, the re- 
mote and nearby zones are bounded by non-degenerate 
roots of p^ (x) and the corresponding physical irrep of the 
fifth-order PAE M.u,r is isomorphic to physical irrep of 
second-order PAE, denoted in the paper as Sf. We use 
here f to distinguish it from collective Dicke index r de- 
scribing the two-level system. We solve the eigenvalue 
problem of the Hamiltonian Eq. (|15|) in terms of simpler 
algebra Sf making use of the isomorphism. The idea is 
straight forward. The generators M±,Mq are realized 
in terms of the generators S± , S3 of Sf according to the 
isomorphism and substituted into the Hamiltonian. The 
interaction part is expanded into perturbation series of 
power of the operator S3 and the series are diagonalized 
by consecutive unitary transformations. 

To begin with we consider the transformation of Mm,t 
to Sf for the case of remote zones. The dimension of a 
remote zone is 2r + 1. Hence, the algebra Sf is charac- 
terized by f = r. The finite dimensional irrep of Sf is 
isomorphic to the corresponding irreducible representa- 
tion of su(2) algebra. The corresponding transformation 
from the generators of Mjy> to the generators of Sf is 
defined as follows (see |34j ) 
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Spectrum {m} of the operator S3 belongs to the interval 
[— r, r] so the argument of the square root function in 
Eq. (|2"5|l has positive valued spectrum in the remote zones 
(M > 2r). Expanding the root function with respect to 
(S3 — h) we obtain the perturbation series with smallness 
parameter a = * 1 . 

It is worth to notice the connection between new op- 
erators S± , S3 and the physical operators of the model. 
From Eas. ljlTH) . (|2l"1) . and IJ2"5|) it follows that in remote 
zones 
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Notice that the subspaces corresponding to the remote 
zones do not contain the vacuum state of the field. It is 
also worth mentioning that the matrix representation of 



the operator 



•y/ata+l 



a is dn^n+i in any remote zone. This 



operator has been considered before in terms of phase 
operator |3ll3ll37l|. 

We turn now to the nearby zones M < 2r. Notice 
that the dimension of nearby zones is q 3 — q\ = M + I , 
and therefore f = 4f. Hence, we obtain the following 
realization 

M = § - S 3 , M_ = (M+)t, 

M + = S_ ^4t_m +1 _g 3 (1 + jk(M + 1 _ §3)) .(26) 

Since all the eigenvalues of the operator S3 belong to the 
interval [— f, r], the argument of the square root function 
does not have zero eigenvalues in the nearby zones. The 
realization of S? through spin and boson variables is then 
given by 

S3 = — a' a, S+ — ; ^^S+ a. (27) 

2 + Vr + 1 - S3 V ' 

Notice that the nearby zones do not contain the eigen- 
vector \r,r) of S3. 



IV. DIAGONALIZATION PROCEDURE 

We start with the representation of the Hamiltonian 
Ea. (|15|) in the form of series of the operator a(S3 — 4), 
where the constant a is a smallness parameter being spec- 
ified below. As we already mentioned, to construct the 
series we expand the square root function in the realiza- 
tions Ea. (|25|l and Eq. l|2()|) depending on the zone under 
consideration. The smallness parameter has the form 



M-r + i : 



4r-M+l ' 



M>2r 



M < 2r 



(28) 



An accuracy necessary to observe all the interesting ef- 
fects described in the paper, is provided by first three 
terms in the expansion of the effective Hamiltonian 
Eq. (|15|) . Up to second order with respect to the smallness 
parameter a it reads 



H « uj p (M + f - r) + A(S 3 +r-r) + 



S++S- 



KO 



S+ + S- 



(**-i)) 



f (5. 



S+ + S- s 3 - 



.(29) 



Here we used a simple relation S3 = S3 + r — r. To prove 
it one should utilize the definition of f in the remote 
or nearby zones given in the previous section, and the 
relations Eas. (|2"5|> and (|2Tjl . The parameters in Eq.J2JJ) 
are defined as follows: 



l + f^(M-r+±), M>2r 



7 = 



k 2 AI+1 
fel 2 ' 



(30) 



M < 2r 



k = ti7^, 



/3i = «(l + lf^)^ 2 =« 2 (|-lf^) 



(31) 



Above, considering the irreducible representations of the 
fifth-order algebra MA/, r we noticed a limitation on the 
value M provided that the algebraic approach for the 
Tavis-Cummings model is applicable. In other words, 
relatively accurate approximation of the structure poly- 
nomial Pb(x) by the parabolic polynomial of Sf in the 
interval between two roots where P5 (x) takes positive val- 
ues is achievable if qo doesn't belong to the interval. Now 
we are able to define the limitation of the algebraic ap- 
proach more rigorously. We find from Eas.l|H() |l .l|Hl |l that 
7 - 0(1) and p n ~ a n if 



M + l < 



,*i. 



(32) 



The Hamiltonian Eq.J35J) may be rearranged into the 
following form 



3 r 



H& = c + as 3 + 

%■ (§ a 3 a + §J 3 ) - fa (s 3 S x S 3 + JS X )] ,(33) 



where S x — —^ — according to the su(2) algebra nota- 
tions, and Co = uj p (M + 4r — r) + A(r — r) is a constant in 
each irreducible representation because M is the Casimir 
operator. 

The first three terms in Eq. I|33|l are linear with respect 
to the generators of Sf and can be diagonalized by well- 
known su(2) unitary transformation Uq = e"' s »' corre- 
sponding to rotation of the quasi-spin vector {S x , S y ,S 3 ) 



about y-axis, where S y 
from the relations 



The angle ^o is found 



cos(V>o) 



n. 



sin (i/>o ) 



O.i 



n R = Va 2 + fc 2 , (34) 



where we introduced a notion of nonlinear quantum 
Rabi frequency Q,r. Hence, after the transformation 
UqH^Uq 1 a zero order contribution into Eq. l|33|l with 
respect to a reads (Co + Q,rS 3 ) that justifies the name for 
the constant £Ir. It describes the frequency of rotation 
of the quasi-spin f of the atom-field quantum system. 

Applying two additional unitary transformations 
U\ , U2 discussed in Appendix [S] we diagonalize the 
Hamiltonian up to second order with respect to a. The 
diagonalized operator takes the form: 



Hfl g {Sz) = {U 2 UiU )HM{U2U x Uo)- 1 = C + tl R S 3 - ^^ (zS 3 2 - r(f + 1)) + 



'*)'£* 



4A -9A z fc^+4fc 



4 s 2 



-A ! ---.A--/,---2fc' - ( - + ±) + 1 Al+A^+fel 



I 2 A- 2 



01 k 2 q 






-Si 



2 -K^f(f- 



1 A A -k A 

2 at 



+ o{a 2 ). 



(35) 



In Fig. 01 we compare the second order solution Eq.JSHJ) 
with the exact numerical diagonalization of H Eq. i|15|) 
and find that the approximation is very accurate and 
doesn't significantly depend on the number M of the par- 
ticles in the condensate. Basically, as higher the ratio 

2r 1S as better the analytical solution Eq. (|35|1 be- 
comes. The first order correction to the spectrum (of 
order (3\ in Eq. (J2SJ)) is proportional to the detuning A 
and vanishes at the point of exact resonance. Therefore, 
the relative error decreases significantly only if the second 
order correction is taken into account. One can under- 
stand from Fig. 0] that at the points where zero and first 
order errors increase to 50% our second order correction 
provides accuracy even higher than 95%. The plot de- 
picts typical behavior of the second order solution in the 
nearby zones. According to the condition Eq. llol'l) the 
approximation diverges when we approach too close to 
the point Ifci/fel — 1 = 100. Therefore, we restrict our 
initial state of the system to the maximum of 60 exci- 
tations. It is justified in the next section from physical 



point of view. 



V. SUBPOISSONIAN DISTRIBUTION IN 
PHOTON STATISTICS 

To study photon statistics we have to construct the 
time evolution of corresponding field operators. It is well- 
known that the simplest characteristic of subpoissonian 
nature of the photon distribution is the Fano-Mandel pa- 
rameter Q(t) defined as follows 



Q(t) 



((ota) 



(ata) 2 



(oto) 



1. 



(36) 



Our unitary transformation approach allows to represent 
the averages of field operators in Eq. 136|) in the following 
convenient form: 



(A(t)) = (<f> \e lHt Ae- tm \$o) 



{QoPo 1 (f&SUoAU^e 



,12) 



0b l*o>, 



(37) 
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FIG. 4: Relative error, 8E{M) - max* 



fi(Af)-E<' 1 )(Af) < 



for 



AE(M) 

three different orders of accuracy (n = 0, 1,2). Here, E(M) 
is the exact spectrum of the Hamiltonian Ea. l|15^ calculated 
numerically in the irrep specified by given M and r values, 
AE(M) is an average splitting between the spectrum levels, 
and E^(M) is the spectrum calculated by Eg. 1351 in n-th 
order with respect to a. For the condensate, we choose N — 
1000, hence r — Af/2. For the field, the parameters A/ui p — 
2.4E-8, kx/ujp = 3.04E-7, k 2 /cj p = -3.01E-9 are taken 
according to the consideration in the previous section. 



where |$o) is an initial state of the atom-field system 
and A is an arbitrary operator. In Eq. I|37[) we left only 
zero order terms with respect to a, i.e. Ui t % — * 1, but the 
diagonalizcd Hamiltonian H^iS^) is kept up to second- 
order terms because the time intervals of interest may be 
relatively large. Applying the relation a) a = M — S3 — f 
it is straight forward to find the time dependence of A = 
a' a from Eq. (|37[l . The result reads 



iH { r> trr f rr-1 -iH\> t 

e d.a S u Q a " aU e dza <> 



M 



_k_ ( 5 +e «[^(5 8+ i)-^(S,)j +hc 



(38) 



Taking the square of the expression Eq. (|38|l one easily 
finds corresponding formulae for A = (a^a) 2 . To cal- 
culate the parameter Q(t) we have to average these op- 
erators over unitary transformed initial state, i.e. over 
Uq\&o)- The transformation Uo is studied in detail in the 
theory of the su(2) algebra. If we succeed to represent 
the initial atom-field state |4>o) in the basis of eigenstates 
|m,f) of the operator S3 then the corresponding ampli- 
tudes -Bm.Ai in the expansion \rh,f)d r = Uo\m,f) — 
J2mi=-f Bm,mi\rh\,f) can be taken from a textbook. 
The state |m, f)d r is usually called generalized atom-field 
dressed state. 

In the paper we choose a typical initial state of the 
system, which seems natural for the experiments with 
sodium BEC. The probe field is prepared in a coherent 
state with average number of photons no — 25. The 



atoms are prepared in completely symmetrized unexcited 
eigenstate of the operator 5 3 , i.e. Jm ,r = 7V/2), with 
zero number of excitations (mo + 4r = 0) . Hence, 



|* > = |-£,£>® E^- na! §\n} 



n=0 



-n ^(L US, _ 



2 ' ' II- 



(39) 



The last equality follows from the relations Eqs. (|21|l and 
(|27|l . because for completely unexcited atoms the total 
number of excitations M is equal to the number of pho- 
tons n. Having represented the initial state in the basis 
of the eigenvectors of the operator S3 we can immedi- 
umtely calculate the Fano-Mandel parameter Q. Below 
in Fig. we plot the evolution of Q(t) and the aver- 
age number of photons for such initial state. Regions 
where the parameter Q(t) takes negative values corre- 
spond to subpoissonian distribution in statistics of pho- 
tons in the probe pulse. First, we notice that the mini- 
mum of Q(t) is reached after the second Rabi oscillation, 
which is described by corresponding Rabi frequency fi^ 
given in Eq. I|34|l . The frequency depends on the number 
of excitations M. Being calculated for the average initial 
number of excitations M = hq it gives the period of Rabi 
oscillation Tr — 2tt/Q,r — 0.12ns, which is confirmed by 
Fig. (see the dynamics of the average number of pho- 
tons). The time interval where the probe pulse obtains 
maximum squeezing in fluctuations of the number of pho- 
tons at, is approximately 0.5ns that is much smaller than 
the relaxation times in BEC (T re i = 16.3ns). This time 
scale satisfies the adiabatic condition imposed at the be- 
ginning. The dispersion of photons in the initial coherent 
state is equal to the average number of photons, i.e. to 
no = 25. It means that the probability to find more than 
2no = 50 photons is negligible and the restriction from 
above on the number of excitations M < l^i/fel ~ 100 
in the system is fulfilled as well. 

We can notice in Fig. [S] that the maximum squeezing 
is relatively large, Q « —0.6. It is observed after two cy- 
cles of almost complete absorption and re-emission of the 
photons by sodium BEC that is represented by the plot 
of the average number of photons in the probe pulse in 
Fig. |3J The quantum effect is provided by the correlation 
in the atomic system, which is transferred to the field af- 
ter two complete Rabi cycles of the collective atom-field 
quasi-spin S3. 

The inline plot of the parameter Q(t) at large time 
scale shows that the effect has very short life-time and 
only appears at the beginning. The regular collapses 
and revivals are observed due to the interference between 
Rabi oscillations with different frequencies. The effect 
is well-known in the one-atom case, and the collapse 
and revival times were estimated in the single-particle 
model|33. We will show now that due to the nonlinear 
interaction between sodium atoms and the probe pulse 
the maximum squeezing of photon statistics is achieved 
much faster than in the regular Dicke model. The fact 




-(1) 

(2) 



t(ns) 



»"f 



10 20 30 



plays crucial role if one takes into account different mech- 
anisms of the decoherence and relaxation in the system. 
For instance, the adiabatic approach developed in the pa- 
per would break down without the nonlinear terms in the 
Hamiltonian. 



t(ns) 



FIG. 5: The Fano-Mandel parameter Q(t) and the average 
number of photons for initially unexcited sodium BEC and a 
coherent probe field with no = 25 average number of photons. 
All the parameters are the same as in Fig. [I] The solid line (1) 
on the plot of Q(t) is calculated in second order of a according 
to Eq. IETF! . The dotted line (2) is calculated by Eqs. JIHt and 

ED. 



Since we are interested in a short time-scale dynamics, 
we can neglect higher order corrections to the spectrum 
of the Hamiltonian Eq. H35[l keeping only zero order terms 
in Eq. $£Q. Then, we find 



(*r 



(oto(*)) = (* |(M - f)(l - cos(^i))sin(^o) 2 |*o) + <$o|a f a (cos(^o) 2 + sin(^o) 2 coa(Sl R t)) |$ ), (40) 

<at a (t) 2 ) = ($o|(a f a) 2 [cos(^o) 4 + cos(V-o) 2 sin(^ ) 2 (3cos(^i) - 1) + \ sin(^ ) 4 (l + 3cos(fi fl i))] |$ ) 

+ ($o|a f a [(M-f)sin(V'o) 2 sin(^-t) 2 (5 + 3cos(2V'o) + 6sin(V'o) 2 cos(fi fl t))] |$ ) (41) 

,| [§ sin(V-o) 2 sin(^t) 2 (3f(f + 1) - (Af - r) 2 + (f(f + 1) - 3(M - f) 2 )(cos(2^o) + 2 cos(O fl i) sin(^o) 2 ))] |$o>- 



Notice, that M, cos(t/jo), sm(tp ), r, Qr are not c-numbers 
but Casimir operators, so they have to be averaged over 
the initial state as well as the operator of the number 
of photons (see the example in Appendix |B|) . The equa- 
tions Eqs. I|40|) and (|41(l are derived in the assumption 
that the initial state of atoms have zero dipole momen- 
tum. But the inversion in the atomic system is allowed to 
have nonzero values. Our initial state satisfies the condi- 
tion and the comparison of second order and zero order 
calculations is provided in Fig.[S]for Q(t). At short time 
scales the agreement is very good. 

In a similar way as it was done by Eberly and co- 
workers [38| , we calculate the characteristic times of col- 
lapses and revivals in the system by the saddle-point ap- 
proach. The details are provided in Appendix [Bj The 
revival time r rev is determined by the difference between 
two adjacent Rabi frequencies Qr with M = no and 
M = no + 1. The collapse time r co ; has two different 
asymptotics. In the case of very small ki 



JX) 



/4(3r-TO + l)fc 2 + A 2 



kfy^ 



(42) 



TABLE I: Characteristic time scales for three different com- 
binations of the coupling constants in sodium BEC. 





ki/u> 


fe/w 


T R (w) 


r£! (ns) 


TJ2 m 


Trev (ns) 


a) 


3E-7 


-3E-9 


0.12 


41.1. 


2.1 


19.9 


b) 


3E-7 


-3E-10 


0.10 


41.1 


20.8 


138 


c) 


3E-6 


-3E-10 


0.01 


4.1 


20.8 


34 



and for relatively large nonlinearity but small 1/V 



.( 2 ) 



'col 



w 



n r 



(43) 



In Fig.EJwe demonstrate a good agreement between these 
zero order estimates for different ratios fe/fci, an d the 
numerical results. It is worth to notice that according 
to Eq. (|43|l the collapse time decreases as k^ 1 allow- 
ing to achieve significant squeezing before the relaxation 
processes take place. As we showed the nonlinear term 
provides additional control on the time-scale of quantum 
effects observed in the system. 
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FIG. 6: The average number of photons for three different 
combinations of the coupling constants fci, &2. All parameters 
are the same as in Fig. [I] The coupling constants and the 
characteristic time scales are given in Table [I] 



VI. CONCLUSION 



close to the frequency of corresponding atomic transi- 
tion. It was demonstrated that the nonlinear refraction 
index can reach extremely large values. The regions of 
its negative values and nonlinear dependence of the ab- 
sorption on the pulse intensity provide a possibility to 
achieve almost complete compensation of the losses and 
the dispersion in the BEC. 

Applying the formalism of polynomially deformed 
su(2) algebras we analyze complex quantum dynamics of 
the probe pulse in a coherent ensemble of sodium atoms. 
It was shown that the nonlinear effects provide a nec- 
essary window in the relaxation mechanisms where the 
probe pulse in the regime of " slow" light exhibits nonclas- 
sical properties in statistics of photons. We demonstrated 
that the poissonian statistics of photons in the coher- 
ent state of the probe pulse can be significantly squeezed 
within short period of time to highly subpoissonian val- 
ues due to collective interaction between coherent atoms 
and the field. These results open new perspectives in 
generation of nonclassical atomic or field states in Bose- 
Einstein condensate controlled by external classical light. 
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APPENDIX A: UNITARY TRANSFORMATIONS 

After the transformation Uq = e l ^° Sy the Hamiltonian 
Eq. (33$ reads 



J 



UoH^Uo 1 = Co + n R S 3 - k^ [cos(2^ ) (§ 3 S X + S X S 3 ) + sin(2^ ) (§$ - S*) 

(cos(3?/>o) + cos(>o) sm(ipo) 2 )S 3 S x §3 - (sin(3'0o) - sin(^ ) cos(tp ) 2 ) S X S 3 S : 
+ sin(V- )S 3 (cos(Vo) 2 (4 2 ~ 1) + l) + cos(i> )S x (sin(Vo) 2 (^ - 1) + ± 



K 2 



(Al) 



To diagonalize the operator Eq. i|Al() in first order of j3\ ~ a we apply second unitary transformation 



U\ = exp 



\ sin(^o) ({S 3 S y + S y S 3 ) cos(2V> ) - \(S X S V + S y S x ) sin(2</> ) 



(A2) 



It is obvious that the transformation U\ applied to the zero order term in Eq. (|A1|I gives the first order term in the 
original Hamiltonian but with an opposite sign providing that it is cancelled. Thus, the Hamiltonian UiUoH^U ~ 1 Ui 1 
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is diagonal up to second order of a. It has the form 



UiUoH^U^U^ 1 = Co + Q R S 3 - fc% sin(^o) cos(^o) (s§l - r(f + 1)) 
(%) 2 sin(^ 



cos(2^o) 2 S 3 (4S$ - 2r(f + 1) + ±) + 22^1! £3(552 _ 5f (f + 1) + ") 



"( 



95*352; S*3 + 5 X 



6f(f+l)-13 1 



+3 sin(2V»o) S X S 3 S X - - 2 ,„„ a „ x „ a , ^ 

(cos(3^ ) + cos(^o) sin('0o) 2 )535 a; 53 - (sin(3?/>o) - sin(^ ) cos(^o) 2 )5z5 3 5 ; 
W>o)5 3 (cos(^ ) 2 (5 3 2 - 1) + i) + cos^ )S x (sin(^o) 2 (5 2 - 1) + ± 



sin 



(A3) 



Notice that we only keep terms of order a or a 2 . Therefore, the last term in Eq. (|A1|) stays intact after the transfor- 
mation. One can see that the first order correction vanishes in the case of exact resonance A = 0. So, it is necessary 
to consider second order terms to take into account significant effects connected with the nonlinear dependence of the 
energy splitting on the value of collective atom- field quasi-spin vector 5 3 . 

To diagonalize the operator Eq. QA3J I up to the order a 2 we apply the third unitary transformation 



lh 



exp 



& sin(^ )) {| sin^o) 2 ^^ + S V S 3 S X ) + 2^|M (§SjJ - 85 3 5,5 3 + f -^~S v - %S V ) } 
2cos(3V>o)5 3 5j / 5 3 + \ cos(ip )S y + \ sin(^o) cos^o) 2 ^.^ + S x S 3 S y ) 
_sjn(^ol { g v§3§x + §J 3 S V ) - cos(V-o) sin(^ ) 2 5\(§5 2 - 2 r(f + 1) + f ) 



— «4r sin(V'o 



.(A4) 



Consecutively applying these three transformations 
C^o,i.2 to the Hamiltonian Eq. i|33|l we obtain the diag- 
onal operator Eq. (|3"5l) . 



APPENDIX B: THE COLLAPSE AND REVIVAL 
CHARACTERISTIC TIMES 



where 



f(x) = x(l — lnx) 
+ihjii + ( 2r + 2 - T) {**£*■ + k ^f- ( B4 ) 



According to the zero-order solution Eq. I|4(J|I a dynam- 
ics of the average number of photons for the initial state 
Eq. (JjSIJ has the form 



" nn^e-^ kj + 2A 2 + fc 2 cos(^T^t) 
n W = 2. —. 2(fc 2 +A 2 ) ' (B1) 

where based on Eqs. (|30() and i|31|) we defined 
n + 1 



k\ + fc 2 



!L ^) V2(4r-n+l). (B2) 



Deriving the expressions we used the fact that our initial 
state belongs to the nearby zones (f = M/2) and the 
atoms are completely unexcited (M = n, mo = — r). For 
simplicity of the analytical expressions below, we repre- 
sent the results taking niQ = —r. But they can easily be 
generalized for arbitrary initial number of excitations in 
atomic subsystem. Since we are interested in the time 
dependence of nit), the constant terms in Eq. (|B1|) can 
be discarded. We estimate the sum of the time depending 
term using the saddle-point method, which has been ap- 
plied to calculate the one-atom model [38| . Denoting the 
time dependent term in Eq. IjBlfl by w(t) we can write 



w(t) 



2vr2(fc 2 +A 2 ) 



Re e 



-n +n f(n/n ) 



dn, (B3) 



Using the approach we find 



w(t) 



n k 



m + a 2 ) 



Re 



,-n Q +n f(xo) 

~7rW 



(B5) 



where the point xq is the saddle-point of the analytical 
function f(x), i.e. /' (xo) = 0. It follows from the defi- 
nition that the saddle-point xq depends on the time. As 
it was explained in detail in the paper [38|. the collapse 
time t co i is roughly estimated by the condition 



|1 - Ref(x )\ = ^/no■ 



(B6) 



For the moment r co / when the condition Eq. I|B6|I is ful- 
filled, the exponent in Eq. I|B5(1 becomes very small and 
the envelope of the sinusoidal oscillations "collapses". 
For t = it is plain to see that Xo = 1. 

Expanding the solution for xq in the vicinity of unity 
and substituting it into Eq. (|B6|) we obtain relatively 
lengthy expression. If we assume that the nonlinear sus- 
ceptibility k<i is the dominating smallncss parameter, and 
expand the result with respect to ki — > 0, we find 
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/4(3r - m + l)kf + A 2 



4fc. 



fc 2 (8r - n + l)(4r + 1) - (2r - n )A 2 



Assuming that the leading smallness parameter is 1/r we obtain an expression for relatively large nonlinearity fc 2 



(B7) 



r (2) - 

' col — 






8/cf + 6fc 2 fc 2 (2n + 1) + 4fc3fc 2 (2n + 3) + fcifcf (f£ + 6n + 1) + fc|(n - ^(n - I)) 

2r/c 2 (2fc 1 +A; 2 ) 3 ' 

I 



(B8) 



Even for fc 2 = in the Eq. (|B7|I our results can't be 
straightforward compared with Eberly's paper because 
we consider the dynamics in the nearby zones, which do 
not exist in the simple one-atom model. But, it is not a 
problem to obtain corresponding results for the remote 
zones, which are more relevant to their work. We also 
want to notice that the solutions Eq. (|B7|I and (|B8(1 pro- 
vide additional information about the dependence of col- 
lapse times on the number of particles in BEC, which is 
completely opposite in these two asymptotics. 

The revival time is estimated by the period when the 
phases of oscillations of neighboring terms in Eq. I|B1|) 
with frequencies k n differ by 2tt. The difference is es- 



timated for the dominant sinusoids with n — no and 
n = no + 1. The revival time reads 



2tt 



(B9) 



lkl o+1 + ^-^/kl a +^ 
The Rabi oscillation period Tr in Tabled] is estimated as 



Tr 



2tt 



V k l, + ^ 
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